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ABSTRACT 

We present Chandra HETGS spectroscopy of the Be/X-ray binary 1A 0535+262 obtained during 
the 2009/2010 giant outburst. These are the first CCD grating spectra of this type of system during 
a giant outburst. Our spectra reveal a number of lines including a narrow Fe Ka emission line with 
a FWHM of ~ 5000 km s _1 . For the first time, we detect the presence of a highly ionized outflow 
in a Be/X-ray binary. Assuming that the line is He- like Fe XXV, fits with a simple Gaussian imply 
an outflow velocity of ~ 1500 km s -1 . However, self-consistent photoionization modeling with xstar 
suggests that Fe XXIII-XXIV must also contribute. In this case, an outflow velocity of ~ 3000 km s _1 
is implied. These results are discussed in the context of the accretion flow in Be-star, neutron star, 
and black hole X-ray binaries. 

Subject headings: accretion, accretion disks — winds — neutron star physics — X-rays: binaries — 
X-rays: individual (1A 0535+262) 



1. INTRODUCTION 

X-ray binaries are divided into a number of differ- 
ent categories based on their observational (e.g. X- 
ray/optical) and physical characteristics . The primary 
division is made based on the mass of the secondary 
star in the binary, with systems for which M2 % 1 Mq 
being classified as low mass X-ray binaries (LMXBs), 
1 ft M2 > 10 M classified as intermediate mass X- 
ray binaries (IMXBs) and those with secondary masses 
greater than 10 Mq classified as high mass X-ray bina- 
ries (HMXBs). Approximately 300 X-ray binary systems 
have been detected in the Galaxy, with more than half 
of these being LMXBs. Of the approximately 115 known 
HMXB syst ems, Be/X-ray binari es comprise the major- 
ity (~ 70%; Ou!iEan[200l[2003). 

Be/X-ray binaries comprise of a neutron star primary 
in a wide eccentric orbit with a Be-star secondary (no 
black hole /Be-star binaries have bee n discovered to date 
- see iBelczvnski fe Ziolkowski|[2009|) . X-ray outbursts in 
these systems typically occur during periastron passage 
of the neutron star close to the circumstellar excretio n 
disc of the Be-star, e.g. lOkazaki fe Negueruelal (|2001j ). 
However, the X-ray outbursts have been observed to oc- 
cur in two distinct flavors, (i) Type-I: Regular outbursts 
which occur during periastron passage, generally repeat- 
ing on the binary orbital period, and (ii) Type-II: Also 
known as giant outbursts likely arising from a dramatic 
expansion of the circumstellar disc. In contrast to type-I 
outbursts, type-II outbursts do not have a strict phase 
dependence and typically last much longer than type-I 
bursts, with the outburst time typically being measured 
in months, e.g. see Fig. [TJ 

Cyclotron absorption lines have been detected in the 
X-ray spectra of a number of Be/X-ray binaries allow- 
ing accurate constraints to be placed on the neutron star 
magnetic field. Typically, magnetic flux densities in the 
range ~ 10 11 — 10 13 G are m e asured, e.g..lNakaiima et all 
l|200fih : iTerada et all l|200"fih : Iwilson et all (|2008f ) . This 
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opens up the possibility of studying the effect of the neu- 
tron star B-field on the accretion flow. In these HMXBs 
the magnetic field is much larger than the B-field of the 
neutron stars in the much older population of neutron 
star LMXBs, e.g., a magnetic flux density of ~ 10 8 G 
has been measured in the accreting m illisecond pulsar 
S AX J 1 808. 4- 3658 IjCackett et aUl2009h . 

iReigl (|2008l ) have investigated the timing properties 
exhibited by sample of 4 Be/X-ray binaries, which 
displayed a type II outburst that was observed by 
RXTE. A number of similarities with the well stud- 
ied behavior of LMXBs were detected when analyzing 
the colour- colour (CD) and hardnes s inten s ity diagrams 
HID) , e.g. lHasinger fe van der K lis (19891): Ivan der Klii 
2006). The main similarities are as follows (i) Similar 
spectral branches in the CD/HID as those observed in 
LMXBs, (ii) smooth motion through the CD/HID, (hi) 
larger variability at low intensities, (iv) power spectra 
that can be described by a small number of Lorentzians, 
and (v) flat topped noise at lower frequencies in the hori- 
zontal branch that turns into power-law noise in the diag- 
onal branch. A number of differences were also observed, 
namely (i) Different patterns in the CD/HID, e.g. a low 
intensity soft state, (ii) slower motion through the var- 
ious patterns (weeks/months as opposed to hrs/days in 
LMXBs), (hi) characteristic timescales that are an order 
of magnitude longer, and (iv) no apparent correlation 
between the power spectral parameters (QPOs, RMS) 
and the mass accretion rate. A number of features also 
hint at a B-field dependence of the timing properties, e.g. 
power-law noise in the HB branch in V0332+53 but flat 
topped noise in 4U01 15+63 and EXO 2030+375 whose 
B-fields are 2-3x lower than that of V0332+53. These 
observations demonstrate that the accretion flow in the 
Be/X-ray binaries is fundamentally similar to that in the 
more thoroughly studied LMXBs. 

Simultaneous observations at optical/NIR and X-ray 
wavelengths have shown substantial changes in the 
structure of both the accretion and circumstellar discs, 
with both discs observed to disappear at times. This 
changing disc str ucture is responsibl e for the varying 
outburst phases ( Wil son et alj I2002D . The observed 
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Fig. 1.— SWIFT BAT hard X-ray lightcurve for the Be/X- 
ray binary 1A 0535+262 from 2005 January to 2010 January. A 
number of regular outbursts (type-I) are detected in addition to 
the giant outbursts (type-II) in 2005 and 2009. The inset shows 
a zoom in on the 2009 giant outburst, the lines mark the times of 
our Chandra HETGS observations (see text). For reference a flux 



of ~ 0.22 counts s 
pulsar. 



cm is normally measured from the Crab 



changes in the accretion and circumstellar discs have 
been successfully explained via the 'truncated disc model' 
( Okaz aki fc Negueruelal I2001D . Here the circumstellar 
disc is truncated due to the tidal interaction with the 
neutron star. Detailed studies of EXO 2030+375 found 
that the circumstellar disc is likely truncated at the 4:1 
tidal resonance radius. This is close to the critical radius 
at periastron. If the disc extends to, or beyond this ra- 
dius, matter may escape through the inner Lagrangian 
poi nt and form an accretion disc a round the neutron star, 
see lOkazaki fe Negueruelal (|2001l) for further details. 

1A 0535+262 is a 103s X-ray pulsar in orbit aro und a 
09.7 - BO Hie companion star ((Steele et al.lll998| ). The 
~ 111 day orbital period is highly eccentric (e ~ 0.47; 
iFinger et al.l[l9 96). A total of 7 giant outbursts have 
been observed from this system to dat e - 1975, 1980, 
1983 1989, 1994, 2005 and 2009, see ICaballero et"aLl 
(2008) for further details. I n response to the 2009 gi- 
ant outburst of 1A 0535 +262 (|Wilson-Hodge et alj|200a 
iRevnolds fc MiHerll2009l ). we requested and were granted 
20ks of Chandra director's discretionary time to obtain 
a HETGS spectrum during this outburst. These are the 
first X-ray grating spectra of a giant outburst from a 
Be/X-ray binary. 

In this paper, we describe observations undertaken 
with Chandra , during the 2009 giant outburst of the 
Be/X-ray binary 1A 0535+262. In §2, we describe the 
observations. We proceed to analyze the data in §3, 
where we find evidence for the presence of a high ve- 
locity highly ionized outflow in this system. In §4, these 
results are discussed and compared to high resolution X- 
ray spectral observations of a number of X-ray binaries, 
and finally our conclusions are presented in §5. 

2. OBSERVATIONS 

2.1. Chandra HETGS 

A giant ou tburst from 1A 0535+262 was detected 
in late 2009 (jWilson-Hodge et all 120091) . the outburst 
peaked at a 15 - 50 keV flux of approximately 6 Crab. 
Unfortunately, as was the case during the previous 2005 



giant outburst, the source was positioned close to the sun 
rendering observations impossible over a large part of the 
outburst. However, the long duration of a typical type-II 
outburst from this system (30 - 60 days) presented an op- 
portunity to obtain the first high resolution spectrum at 
X-ray wavelengths of a Be/X-ray binary in outburst. As 
such, we requested, and were granted, a Chandra DDT 
observation. 1A 0535+262 was observed on 2 occasions 
upon exit from a sun constrained sky position in late 
December 2009. 

In the first observation (obsid: 12066 - 2009.12.28), 
the exposure time was ~ 9 ks, while the exposure time 
for the second observation (obsid: 12067 - 2009.12.31) 
was ~ 11 ks. The high energy transmission grating spec- 
trometer (HETGS) was used to disperse the incoming 
photons onto the ACIS-S CCD array. As the source flux 
was high, observations were carried out in continuous 
clocking mode, with a nominal frame time of ~ 2.8 ms. 
In this data mode the fast readout times are achieved 
by collapsing the data into a single dimension, preserv- 
ing spectral and timing information at the price of losing 
one of the spatial dimensions. 

Spectra and lightcurves were extracted from the Chan- 
dra event lists using CIAO v4.2. Response files were gen- 
erated using the MKRMF & fullgarf tasks and first 
order spectra from each grating were combined with the 
ADD_grating_ORDERS task resulting in a single HEG & 
MEG spectrum for each observation. The spectra were 
binned using the GRPPHA tool to ensure 10 counts per 
spectral bin before exporting to xspec for analysis. 

2.2. RXTE 

1A 0535+262 was also observed by RXTE during thi s 
outburst on an almost daily basis (ICaballero et alf 2009) . 
Data was obtained with both the PCA and HEXTE. The 
RXTE observations were quasi-simultaneous with the 
Chandra observations with pointed observations occur- 
ring within a 12hr period of our observations on each oc- 
casion (Obsid: 94323-05-03-03,-06; PI: Cabellero). Spec- 
tra were extracted in the standard manner using the rele- 
vant FTOOLS. For the PCA, spectra and background files 
were extracted from layer-1 of PCU2 alone, giving useful 
data in the spectral range 3-20 keV. HEXTE spec- 
tra were extracted from cluster-A, while backgrounds 
were extracted from cluster-B and converted to cluster-A 
background files using the FTOOL HEXTEBACKEST 2 giv- 
ing useful data in the 15 - 100 keV band. As a number 
of background features were apparent at higher energies, 
data above 60 keV was ignored in all further analysis. 
Hence, the RXTE spectra resulted in a broadband de- 
tection of 1A 0535+262 in the spectral range from 3 - 
60 keV. 

3. ANALYSIS & RESULTS 

The HEG spectra immediately reveal the presence of 
an emission feature consistent with iron Ka at approxi- 
mately 6.4 keV in addition to an absorption feature at ~ 
6.7 keV. In contrast, outside of a number of line like resid- 
uals in the Si/S region (2.1 - 2.4 keV), which are likely 
of an instrumental nature, there are no obvious spectral 
features in the MEG spectra. Before we examine in de- 
tail the emission & absorption lines in the spectrum, we 

2 see http://hcasarc.gsfc.nasa.gov/docs/xte/xhp_new.html 
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Fig. 2. — Broadband spectra of 1A 0535+262 as measured 
by Chandra and RXTE. The spectra (0.9 - 60 keV) include data 
from MEG, HEG, PCA and HEXTE. The top spectrum corre- 
sponds to observation 1 (red) while the bottom spectrum corre- 
sponds to observation 2 (black). The best fit model corresponds 
to pha* (cyclabs* (diskbb+bb+cutof f pi) ) , sec Tablc[l]for the best 
fit parameters. 

first model the spectral continuum. To this end, we uti- 
lize the RXTE data in combination with the Chandra 
data (HEG & MEG) to constrain the continuum. All 
data analysis takes place within xspec 12.5.0. 

3.1. Continuum model 

In Fig. [5J we plot the broadband Chandra and RXTE 
spectrum of 1A 0535+262. The spectrum was fit with 
a model consisting of a disk blackbody and a blackbody 
plus a powerlaw with a spectral cut-off containing an 
Cyclotron absorption line modified by interstellar ab- 
sorption (pha* (cyclabs* (diskbb+bb+cutof f pi) ) ), as 
this was consistent with the observed spectrum in the 
2005 normal outburst (jCaballero et all l2008t ). The 
column density was held fixed at a value of 4.5 x 
10 21 cm -2 cons istent with the meas ured value towards 
1A 0535+262 dKalberla et alj 120051) a nd with optical 
measurements, e.g. ISteele et al.l (|1998[ ). Allowing the 
column density to vary results in a best fit value con- 
sistent with this value. At high energies a broad ab- 
sorption feature is observed consistent with the ex- 
pected position of the fundamental Cyclotron resonance 
absorption line in this system. The measured ener- 
gies are consistent with measurements from the previ- 
ous 2 005 outburst (jTerada et al.1 120061: ICaballero etall 
I2008D. and that meas ured earlier in the 2009 outburst 
(|Caballero et al.ll2009| ) implying a neutron star magnetic 
field of ~ 4 x 10 12 G. 

The resulting best fit model parameters arc displayed 
in Table [TJ At higher energies, no significant change is 
observed in the position of the Cyclotron lines; however, 
as the background dominates above 60 keV, our limited 
ability to constrain the continuum above the Cyclotron 
line will affect our ability to quantify any small changes 
that may be present. The simple cut-off power-law model 
that we have used to describe the high energy emission 
provides a good description of the data. If this emission 
arises through Compton up-scattering, then the spectral 
break may signal the electron temperature of the hot 
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Fig. 3.— Chandra HEG spectra of 1A 0535+262 taken at the 
times indicated in Fig. [T] Emission lines from Fe Kc? & K/3 are 
clearly detected in both observations and are consistent with no 
change within the errors. In contrast absorption consistent with 
Fe XXV is only detected in the first observation. The exposure 
time was ~ 10 ks in each observation. Note: The data have been 
rebinned for visual clarity. 

corona. The low energy spectral break we measure (~ 
20 keV) is consistent with a low electron temperature 
if the Comptonizing corona is somewhat optically thick, 
i.e. kT c - 10 keV,r > 1. 

Assuming a distance to 1A 0535+262 of 2 kpc and a 
canonical neutron star mass of 1.4 M©, we measure unab- 
sorbed 0.9 - 60 keV luminosities of L Xl = (1.59 + 0.01) x 
10 37 erg s" 1 and L Xa = (1.17+0.01) x 10 37 erg s" 1 , which 
correspond to Eddington scaled luminosities of ~ 9% and 
6.6% for the first and second observation respectively. 

3.2. Emission & absorption lines 

Having characterized the continuum, the spectra were 
then inspected for emission and absorption features. We 
concentrate on the HEG spectra as there are no obvi- 
ous spectral features in the MEG spectra. A number of 
emission and absorption features are clearly detected in 
the iron K region of the HEG spectrum. In particular, 
we detect emission consistent with Fe Ka and Fe K/3, 
while a single absorption line consistent with the pres- 
ence of highly ionized Fe XXV is also detected in the 
first observation alone, see Fig. [3] 

Basic line parameters were obtained by fitting Gaus- 
sians to the detected lines, the parameters of which we 
display in Table [2] In Fig. [3j we display the HEG spec- 
tra for each epoch with the best fit continuum model 
from §3.11 in addition to a number of Gaussian lines. 
The iron K-shell lines are detected in each observation 
and have the same equivalent widths (EW) within the 
errors. The Ka line is resolved and hence we measure 
a line width of ~ 5000 km s _1 in the first observation 
while the measured line width is lower in observation 
2 at ~ 4100 km s _1 , it is nonetheless consistent with 
the first observation within the error bars (see Table [2]) . 
The observed Ka line widths are consistent with velocity 
broadening in, for example, an accretion disk. The mea- 
sured line centroid (6.39+0.01 keV) supports this, being 
inconsistent with the presence of higher ionization states, 
that would broaden the line, contributing to the observed 
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Fig. 4. — Main panel: Zoom in on the 5-8 keV region of the first Chandra observation with the best fit continuum model. The 
emission lines have been fit with Gaussians. A single absorption line remains in the residuals consistent with an origin from highly 
ionized Fe XXV (bottom panel). An XSTAR absorption component was then added to the model, the resulting best fit model gives 
N = 2.0^Q g X 10 21 cm -2 , £ = 4»0jlo 2 er S cm s_ ant ^ v °ut = 3150 ± 600 km s _1 (lcr errors). The absorption line is found to require 
a contribution from Fe XXIII - XXV, an additional line due to Hydrogen like Fe XXVI is also predicted. Inset: Best fit model of the 
absorption lines in 1A 0535+262, where the y-axis is in units of photons cm -2 s _1 keV - 1 , i.e., the predicted incident spectrum. The 
normalization of the Fe K«, K/3 lines have been set to zero in order to emphasize the predicted absorption lines. Note: The data in the 
main panel has not been divided by the detector area, in order to illustrate the absence of instrumental features in the Fe K region, and 
that the data have been rebinned for visual clarity. 



Fe Ka line, i.e. > Fe X. The measured line fluxes (see 
Table© return a Fe K/3/Ka flux ratio of 0.21±0.11 in 
agreement with th e theoretically expe cted value within 
the errors (~ 0.13. iPalmeri et al.ll2003ft . 

An absorption line, at an energy consistent with highly 
ionized Fe XXV is detected in the first observation alone. 
We do not observe any absorption from Fe XXVI in ei- 
ther epoch, nor is there evidence for any additional ab- 
sorption features in either observation. Assuming the de- 
tected absorption line is due to highly ionized Fe XXV, 
we measure a significant blue-shift consistent with the 
presence of an outflowing wind in the 1A 0535+262 sys- 
tem. The measured line position implies a velocity for 
the outflowing material of v out = 1500 ± 1000 km s _1 . 

3.3. Ionized absorption 

To provide a more physical description of the highly 
ionized absorption in the spectrum of 1A 0535+262, 
we generated a grid of models using xstar version 
2.1kn9, and the "xstar2xspec" facility. For an input spec- 
tral form, we used the (unabsorbed) best-fit continuum 
model, extended to cover the 0.1-100.0 keV range. A 
covering factor of 0.2 was assumed, cons istent with other 
X-r ay binaries (e.g iMiller et al.l [20081 on GRO J1655- 
40, IUeda et"aT]|2009l on GRS 1915+105), and consistent 
with the absence of any strong ionized emission lines. A 



maximum density of n = 10 12 cm" 3 was assumed, and 
initially solar abundances were assumed for all elements. 
This xspec table model was then added to the contin- 
uum model in ^3.11 as an additional absorption compo- 
nent, i.e. pha* (mtablejxstar . table}* (continuum) ) . 
The free parameters of this additional absorber are the 
column density, N, the log of the ionization parameter, £ 
and the redshift/blueshift, Vj n / out . 

Fits to the data with this grid of absorption models 
find the following best-fit parameters: N = 2.9^'q x 

10 21 cm" 2 , f = 3.4l° j erg cm s" 1 , and v out = 5000 ± 
600 km s" 1 (uncertainties are la errors). However, this 
model predicts H-like Mg and Si absorption lines that 
are not observed in the data. Moreover, it appears that 
the absorption at ~ 6.73 keV also contains a contribution 
from Fe XXIII and Fe XXIV that is comparable to that 
from the Fe XXV line. 

A model where iron is over-abundant by a factor of 
two provides an improved fit to the data, and does not 
predict observable Mg and Si lines. With such a model, 
the following absorption parameters are measured: N = 

2.0to:8 x 1021 cm ~ 2 ' £ = 4 - ±a2 er S cm s_1 i and v out = 
3150 ± 600 km s -1 (uncertainties are la errors). This 
model provides a fit both to the entire spectrum and to 
the ionized Fe K region that is more reasonable. In the 



Chandra grating spectroscopy of 1A 0535+262 



5 



Fig. [4] inset, we display the best fit model containing the 
xstar absorption component but with the normalization 
of the emission lines set to zero for clarity. The incident 
model spectrum is plotted to illustrate the contribution 
from the individual Fe ions. The absorption line is found 
to be the sum of absorption from a number of highly 
ionized Fe ions i.e. Fe XXIII, Fe XXVI, Fe XXV, see 
iKallman etldl (|2004l ) for details. 

We emphasise, that it is impossible to produce absorp- 
tion solely due to Fe XXV. In each model where ab- 
sorption due to Fe XXV is produced, we also obtain a 
contribution from the two preceding Fe ionization stages 
(XXIII, XXIV). We can estimate the column density in 
each ion via the standard relation bet ween the line equiv- 
alent width and the curve of growth (|Kotani et al.ll2000L 
2006). Our best fit model (see Fig. 2|) implies column 
densities of - 10 17 cm" 2 for Fe XXV and Fe XXIII, while 
the Fe XXIV column density is less than a third of this. 

The best fit XSTAR model also predicts absorption due 
to H-like iron. This absorption line is not detected in the 
data due to a combination of the narrow line width, low 
S /N and relatively low spectral resolution at this energy 
(~ 7.04 keV). This line will contribute to the red wing of 
the Fe K/3 line, resulting in a larger line EW than mea- 
sured in ^3.21 In turn this will effect the measured Fe 
K/3/Ka flux ratio, increasing it from the value measured 
earlier but nonetheless remaining consistent with the ex- 
pected theore t ical v alue within the errors, i.e., ~ 0.13 
iPalmeri et al.l (|2003| ). However, further observations are 
required to test for the presence of this Hydrogen like 
absorption line. 

4. DISCUSSION 

We present the first high resolution X-ray spectrum 
of a type-II outburst from the Be/X-ray binary 1A 
0535+262. In the spectra we detect emission lines from 
Fe Ka & K/3 in two separate observations separated by 
3 days, which do not display any significant evidence for 
variability between the observations. In the first epochs 
observation alone, we find a blue-shifted absorption line 
consistent with an origin in highly ionized Fe. Detailed 
modelling implies an outflow velocity of approximately 
3000 km s _1 . We now proceed to discuss the observed 
spectrum in more detail. 

4.1. The emission line spectrum 

We can use the observed lines to place constraints on 
the accretion flow geometry. Firstly, we consider the 
observed fluorescence Fe Ka emission line. Assuming 
Keplerian rotation, the line FWHM implies an emis- 
sion radius of r = 7200^35™ km in observation 1 and 

r = 10800ig3oo° km in tn e second observation. This 
emission radius (i.e., ~ 10 4 km) is consistent with the 
Alfven radius given the known magnetic field in 1A 
0535+262 of - 4 x 10 12 G, i.e. r A » 3.7 x 10 3 M~ 2 / 7 km, 
where the mass accretion rate is measured in units of 
~ 10~ 8 M Q yr _1 . Hence, if the Fe fluorescence line orig- 
inates at the inner edge of the accretion disc, which is 
truncated close to the Alfven radius, this implies a mass 
accretion rate < 10~ 8 M Q yr _1 . This accretion rate is 
consistent with the observed X-ray luminosity during our 
observation, Lx ~ 0.09 L<Edd- 

We must also ask if the observed fluorescent line 
strength is consistent with theoretical expectations for 



fluorescent line formation. iGeorge fe Fabian] (|1991l) pre- 
dict the fluorescent line eq uivalent width to be EW ~ 
(£l/2n) * 180 eV. Following IChen fe Halpernl (fl989h . we 
find the solid angle subtended by the disc to be f2 < 0.25. 
The measured EW for the 6.4 keV iron line is ~ 25 eV 
(see Table [2]), implying a solid angle of ~ 0.15. For com- 
pleteness, we can also ask if the observed line could be 
excited on the surface of the mass donor 09.7IIIe star; 
however, for the known orbital separation at perias tron 
(~ 1 x 10 13 cm, iCoe et al J 12005 iHaieh et alJl200l the 
secondary star will subtend a solid angle over an order of 
magnitude less than that of the disc. Similarly the ob- 
served line width is difficult to produce in a stellar wind, 
for example the typical terminal velocity for the wind 
from a B-star is < 2000 km s _1 , this is inconsistent with 
the measured spectrum. 

Iron fluorescence emission lines have been de- 
tected in num e rous other neutron star X-ray binaries. 
iTorrejon et al.l (|2010( l have recently completed a system- 
atic analysis of Chandra HETGS observations of 41 X- 
ray binary systems (10 HMXBs & 31 LMXBs). They 
find a narrow Fe Ka emission line to be present in all of 
the HMXBs studied while it is only present in a small 
fraction of the LMXB sample (< 10%). However, the 
width of this line is measured to be low, with a typical 
velocity width of < 800 km s _1 and is found to be con- 
sistent with a wind like origin. It is clear that the line 
width we measure in 1A 0535+262 (v ~ 5000 km s _1 ) 
is not consistent with having the same origin as the nar- 
row Ka lines detected by Chandra in the majority of 
HMXBs observed to date. Nonetheless, we note that the 
line and continuum flux we measure for the Fe Ka line 
ar e consistent with the re lation Fu„ e oc F° Q 7 * as measured 
by iTorreion et al.l ()2010l ) for the sample of systems they 
studied. This is to be expected given the small radius 
implied by the line width (~ 10 km). The equivalent 
width of the Fe Ka line is also consistent with the low 
significance X-ray Baldwin effect' observed in the other 
HMXBs. 

The line we observe in 1A 0535+262 is produced in 
a similar manner to the relativistically broadened lines 
that are observed in many of the neutron star LMXBs, 
i.e., fluorescent X-ray emission from the surface of the 
accretion disc. However, in 1A 0535+262 the observed 
Ka line does not display any evidence of the characteris- 
tic relativistic red wing as expected if the line is formed 
deep in the g r avitat ional well of the neutron star, e.g., 
ICackett et all (|2010h . Instead the line we measure is 
narrow, consistent with a much larger emission radius 
(~ 10 4 km, see above). 

The emission line spectrum that we observe in the 
1A 0535+262 spectrum is fundamentally different to 
that observ ed in the high mass X -ray binary pulsars 
GX 301-2 (|Watanabe etail l2M$\ and 4U 1908+075 
(jTorrejon et alJl2010t ). IWatanabe etaH (|2003h observed 
GX 301-2 for 40 ks with Chandra HETGS and detected 
a fluorescence Ka line with an accompanying Compton 
shoulder. This line is formed in the dense wind from 
the B2 supergiant secondary star. Detailed XSTAR mod- 
elling of this X-ray spectrum found column densities of 
Nh ~ 10 24 cm~ 2 . Analysis of a Chandra spectrum of 4U 
1908 +075, reveals a simil ar Compton broadened Fe Ka 
line (|Torreion et al.|[2010t ). In contrast the iron absorp- 
tion lines detected in 1A 0535+262 are consistent with a 
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TABLE 1 

Broadband continuum fit parameters 



kT diakbb 


norm diskbb 


kT bb 


norm bb 


r 


E C ut 


norm cutoffpl 


Ecyc 


D C y C 


Wcyc 


[ keV ] 


xlO 4 


[ kcV ] 


xl0~ 3 




[ kcV ] 




[ keV ] 




[ keV ] 


0.15+0.04 


5±1 


1.07±0.04 


7.7 ±0.05 


0.54±0.01 


23.0 ±0.2 


0.28±0.01 


42.85±0.8 


0.75±0.05 


14.9±1.3 


0.17±0.01 


3.1 ±0.3 


1.03 ±0.2 


6.3 ±0.03 


0.44±0.1 


17.6±0.3 


0.183 ±0.004 


44.0±0.5 


0.57±0.06 


8.0±1.0 



Note. — Best model parameters for the 1A 0535+262 continuum in the spectral range 0.9 - 60 keV, see Fig. [2] where the best fit model 
is pha* (cyclabs* (diskbb+bb+cutof f pi) ) . All errors are quoted at the 90% confidence level. 



much lower column density in the wind, N ~ 10 21 cm 2 . 

4.2. Wind origin? 

We have measured an absorption component consis- 
tent with a significant outflow in this system. Simple 
phenomenological modelling favours a outflow velocity of 
~ 1500 km s _1 consistent with absorption due to highly 
ionized He-like iron, see £13.21 Detailed modelling with 
the XSTAR plasma photoionization code reveals the ab- 
sorbing component to consist of a number of highly ion- 
ized Fe ions from XXIII - XXV at a common blue-shift 
of ~ 3000 km s _1 , see £13.31 This is the largest wind 
velocity measured from a Galactic X-ray binary to date. 

We can estimate the maximum radius at which the ab- 
sorption occurs via the measured ionization and column 
density from the xstar model, i.e., assuming N=nr, we 
have r = L x /N£, which implies a radius of ~ 10 12 cm. 
This radius is less than the estimated Roche lobe ra- 
dius of t he neutron star close to periastron (Rl ~3.5x 
10 12 cm: iCoe et all [20061: iHaTgh et al.ll200l and implies 
a disc origin for the observed absorption component. 

There are three well studied mechanisms for launching 
a disc wind (i) radiative - radiation pressure may drive a 
moderately ionized wind (jProga et al.ll2008h (ii) thermal 
- Compton heating of the accretion disc by the central 
X-ray source may raise the gas temperature to the es- 
cape velocity and d rive a thermal wind from the system 
(|Woods et alJ ll996). and (iii) magnetic - could be pow- 
ered either by internal m agnetic visco sity (MRI) giving 
rise to a slow dense wind (|Progall2003fl or rigid magnetic 
fields which would allow some material to be acceler- 
ated a long magnet ic field lines giving rise to a clumpy 
wind ( Spruit 1996). Radiative and thermal winds are 
expected to be launched from significantly larger radii 
(e.g., 10 10 cm vs 10 9 cm for a 10M Q BH) and have a 
lower density than a magnetic wind, a prediction that 
may be tested via detailed modelling of the disc ab- 
sorption line spectrum as demonstrated by iMiller et al.l 
2006al 120081 The radius implied by our modeling favours 
a radiative/thermal origin for the large outflow velocity 
measured in 1A 0535+262. 

4.3. Outflows in other X-ray binaries 

Chandra o bservations of GRO J1655-40 ((Miller et al.l 
l2006al I2008T ) have provided the most revealing insights 
into the structure of this accretion disc wind. Detailed 
studies of the wind absorption line spectrum favour a 
magnetic origin with an inner launching radius ~ 10 9 cm. 
The mass flux in the wind was constrained to be at least 
5% - 10% of the total accreted mass and could eas- 
ily be comparable with the mass accreted through the 
disc. Observations of th e stellar mass blac k holes HI 743- 
322 (IMiller et alJ (2006b) and 4U 1630-472 dKubota et al.1 



2007; have also displayed evidence for dense winds 
launched from the inner disc region. Recent observations 
of the Galactic microquasar GRS 1915+105 also display 
blue-shifted absorption features cons istent with an out - 
flow velocity of 150 - 500 km s" 1 (jUeda et al.ll2009ft . 
While observations of the Galactic high mass black bi- 
nary Cyg X-l al so display evidence for a low velocity 
outflowing wind (|Miller et al.ll2005lh The low outflow 
velocity here is interesting as in this system there is a 
large scale wind, due to the O-type secondary star, with 
an expected terminal wind velocity of Voo < 2100 km s . 
In all of the above cases the measured outflow velocity is 
low < 500 km s _1 . 

Similarly, Chandra observations of the neutron star 
LMXB GX 13+1 also measu re a low velocity, highly ion- 
ized wind (|Ueda et al J 12004). Detailed modelling of the 
observed absorption lines revealed the wind to be consis- 
tent with being radiatively driven. Of the sample of X- 
ray binaries observed with the HETGS, only a single sys- 
tem has displayed absorption lines wit h a velocity com- 
parab le to that which we observ e here. iBrandt fe Schulzl 
(|200Cfl ; ISchulz fe Brandt! (120021) presented spect ra of the 
neutron star binary Cir X-l (|Jonker et al.|[2007T) . The Fe 
XXV absorption line was found to exhibit a P-Cygni pro- 
file consistent with an outflow velocity of ~ 2000 km s _1 . 
In this case it is likely that both thermal and radiation 
pressure play a significant role in accelerating th e wind 
to the measured velocity (|Schulz fe Brandtll2002D. In a 
subseq uent Chandra observation of Cir X-l. lSchulz et al.l 
(2008) unambiguously detect a number of absorption 
lines which they identify as highly ionized Fe XXII - 
XXV. The detection of these highly ionized iron lines in 
Cir X-l would support the results of our XSTAR modelling 
which suggest a contribution from multiple Fe ionization 
levels to the measured absorption line in 1A 0535+262. 

The fact that both of the X-ray binaries with the 
largest observed disc wind velocity (1A 0535+262, this 
paper & Cir X-l. ISchulz fe Brand t 2002) contain a neu- 
tron star primary suggests that the additional hard X- 
ray emission originating from the neutron star surface 
and the disc/neutron star boundary layer is crucial to 
accelerate a wind from the disc to the observed velocities 
(v ut ~ 3000 km s _1 ). Deeper future observations will 
no doubt shed light on precise nature of the observed 
absorption lines and the wind acceleration mechanism 
in Be/X-ray binary systems. Studies of Galactic OB 
associations often find that Fe is under- abundant (e.g. 
iDalfon et ail 1200 ID . The fact that our data and models 
would appear to require an over-abundance of iron can 
be tested via future HETGS observations. 

5. CONCLUSIONS 

We have presented high resolution X-ray grating spec- 
troscopy of the Be/X-ray binary 1A 0535+262. For the 
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TABLE 2 

HEG SPECTRAL FEATURES 



Line 


Energy 
[ keV ] 


EW 
[eV] 


Observation 1 
Flux 

[ photons s" 1 cm" 2 ] 


(7 

[ km s" 1 ] 


Energy 
[ keV ] 


EW 

[eV] 


Observation 2 
Flux 

[ photons s" 1 cm" 2 ] 


a 

[ km s" 1 ] 


Fe I Ka 
Fe I K/3 
Fe XXV Ka 


6.39±0.01 
7.061±0.015 
6.730±0.015 


27 ±f 

7+13 
'-3 

10+_? 


(2At° i) x 
(0.5 ±0.2) x 10" 3 
(0.8 ±0.3) x 10-3 


5090+ 1150 
ou»u_ 1040 


6.39±0.01 
7.058 
6.7 


29 Ij 5 
o+16 
S -6 

< 3 


(1.6 ±0.3) x 10" 3 
(0.4 ±0.2) x 10" 3 
< 0.2 x 10" 3 


4140+ 1240 



Note. — Atomic lines detected in the 1A 0535±262 HEG spectrum. All errors are quoted at the lcr confidence level. Note: where lower 
limits are not given, the lower limit is consistent with zero. Likewise parameters without errors were held fixed. 



first time we detect absorption consistent with the pres- 
ence of a highly ionized outflow in this class of X-ray 
binary. When considered together with the previously 
detected X-ray ionized winds in a variety of other X-ray 
binary classes including: black hole X-ray binaries (e.g., 
high-mass - Cyg X-l, low-mass - GRO J1655-40) and 
neutron star low-mass X-ray binaries (e.g. Cir X-l, GX 
13+1) this would support the emerging picture of a con- 
tinued large scale outflow at large accretion rates but in 
a non-jet form. 
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